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Abstract 

Several  systematic  effects  that  can  influence  SBAS  and  GPS  time  transfers  are  discussed. 
These  include:  ionospheric  delays,  multipath,  and  signal-to-noise  ratio.  Estimates  for  the 
magnitude  of  these  effects  are  given. 


MOTIVATION 

Timing  Laboratories  should  pay  more  attention  to  Multipath  (MP) 

•  MP  is  a  systematic  effect  that  can  be  mitigated 

•  Calibration  of  receivers  through  transportation  of  a  unit  does  not  remove  effects  of  MP 

•  MP  is  not  constant  (can  change  with  time) 

•  Code-phase  MP  error  is  not  zero  mean,  hence,  it  can  not  be  eliminated  through  averaging. 


INTRODUCTION 

•  WAAS  has  taken  great  steps  to  reduce  effects  of  MP  (timing  community  can  learn  from  them  [1]). 

•  More  modem  receivers  (usually  dual-frequency  receivers)  are  mitigating  MP.  Most  timing  receivers 
are  single-frequency. 

•  Several  techniques  have  been  proposed  for  showing  the  existence  of  MP  and  mitigating  it. 


WHAT  IS  MULTIPATH  ? 

•  MP  is  the  corruption  of  the  received  signal  by 
one  or  more  slightly  delayed  reflections  of  the 
direct  GPS  signal. 

•  The  effect  within  the  receiver  is  a  distortion 

of  the  correlation  function  and,  therefore,  of  the 
phase,  code,  and  SNR  measurements  [2], 
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RECOGNITION  OF  MP 


Figure  9.  ^(Ll)  Signatures  as  Evidence  of  Multipath  for  PRN  31  at  Naha 

(JO  1J  S4(Z,1)  scale  between  dates) 


(Figure  is  from  [3].) 


CHARACTERIZATION  OF  MP 

•  1  s  pseudorange  errors  are  approximately 
1  meter  at  low  elevation 

•  1  s  pseudorange  errors  approach  0.2  meters 
at  high  elevations 

•  Consistent  with  Multipath  Estimating  Delay 
Lock  Loop  (MEDDL)  receivers 

•  For  timing  accuracy  to  approach  1  ns,  we 
need  to  do  better. 

(Figure  is  from  [1].) 
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Figure  2.  LI  Pseudorange  Errors  for  WAAS  Network 
(5  June  2000,  CSMOOTH  100,100) 
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ANOTHER  CHARACTERIZATION  OF  MP 


Then,  the  total  S4,  including  the  effects  of 
ambient  noise,  is  defined  as  follows: 


where  (•)  represents  the  expected  (or  average)  value  over 


the  interval  of  interest  (60  seconds). 


Figure  3.  Separating  Multipath  and  Amplitude 

Scintillation 
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(Figure  is  from  [4];  CCDiv  =  Code  Carrier  Divergence.) 


TECHNIQUES  TO  MITIGATE  MP 

•  Siting  -  most  critical 

•  Antenna  -  ground  plane,  choke  ring,  adaptive  antenna  array 

•  Delay  Lock  Loop  (DLL)  -  improved  receivers 

•  Signaling  technique  -  longer  C/A  codes  (L5) 

•  Data  processing  -  code  carrier  smoothing,  SNR  measurements,  repeatability,  multiple  receivers 

•  Carrier  smoothing  techniques 

•  Spectral  decomposition 

•  CNMP  algorithm  of  WAAS. 
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HATCH  FILTER 

•  The  Hatch  filter  is  one  of  the  most  well-known  and  simplest  schemes. 

•  It  estimates  the  carrier  phase  bias,  Xc,  by  averaging  the  difference  of  code  and  carrier  measurements: 

nK(n)  =  -  i  [<pc  (/)  -  Tc  (i) ] 

n  ;=1 


where  n  is  the  length  of  the  data  set. 

•  Using  the  carrier  phase  bias  as  given  above,  it  can  be  shown  [5]  for  the  single  frequency  case  that  the 
pseudorange  can  be  estimated  recursively  by: 

Pin)  =  [p(n-l)  +  A <pL1(n)\  +  ±-Tu(n) 


where  p  is  the  smoothed  pseudorange.  This  smoothing  scheme  is  actually  a  moving  average  window  of 
width  k. 


SPECTRAL  ANALYSIS  AND  MP  DECOMPOSITION 


•  Spectral  analysis  of  MP  is  the  fundamental  concept  for  this  technique  [5] 

•  Time  constants  of  MP  errors  can  have  values  from  1  second  to  hours 

•  Effective  decomposition  of  multipath  error  in  the  frequency  domain  is  critical  in  MP  mitigation 
techniques: 

High  Frequency  (HF)  component  (1  sec  <  x<  1  min);  noise  and  fast  fading  due  to  diffuse  multipath 
and  receiver  tracking  errors 

Medium  Frequency  (MF)  component  ( 1  min  <  x  <  8  min)  due  to  specular-dominated  MP  errors 
Fow  Frequency  (FF)  component  (x  >  8  min)  due  to  long-term  bias-like  MP  and  drift  errors. 


WAAS  CODE  NOISE  MULTIPATH  (CNMP)  TECHNIQUE  [1] 


•  Corrects  for  observed  multipath  and  provides  a  noise  estimate  for  the  residual  error  in  the  MP  corrected 
pseudorange  measurements 

•  The  CNMP  algorithm  has  three  major  components: 


Mean  filter 
Mean  error  function 
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Cycle  slip  detection 

•  The  MP  correction  is  computed  by  differencing  the  mean  filter  estimate  from  the  current  multipath 
estimate 

•  The  MP  corrected  pseudorange  measurements  are  carrier-smoothed  with  a  period  equivalent  to  the 
mean  filter  time  constant. 


WAAS  CODE  NOISE  MULTIPATH  (CNMP)  EQUATIONS 

MPL1[t]  =  Pl iW  -  {(1  +  2k)  tpLl  [t]  -  2 ktpL2[t]} 
MPL2[t ]  =  pL1\t]  -  {{2k  +  2)(pLl[t]  -  (2k  + 1)  cpL2  [/]} 


MP  correction: 


MLj  =  MPLj\t\  - MPLj\t ] 


where: 

p  =  pseudorange 
(p  =  carrier  range 
k_  120 2 
~~  (1542  -1202) 
j  =  LI  or  L2 

t  =  time  of  continuous  LI  and  L2  tracking. 

CONCLUSIONS 

•  Care  must  be  taken  by  the  timekeeping  community  if  you  want  nanosecond  accuracy  from  GPS  time 
transfers 

•  Cannot  obtain  nanosecond  accuracy  unless  one  accounts  for  MP 

•  There  are  several  ways  to  mitigate  MP 

•  One  must  carefully  investigate  one’s  locale  before  choosing  a  method 

•  Calibration  techniques  using  a  reference  receiver  should  carefully  investigate  MP  at  the  sites  to  be 
visited 

•  Observations  should  be  made  with  a  dual-frequency  receiver. 


DISCUSSION 


•  Can  one  really  expect  better  than  a  nanosecond  from  GPS  time  transfers? 

•  For  1  nanosecond,  should  other  techniques  be  emphasized? 
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Questions  and  Answers 

TOM  CLARK  (Syntonics):  A  couple  of  comments.  One,  you  mentioned  the  dual  frequency  problem. 
Sometimes  dual  frequency  can  actually  hurt  you  more.  The  normal  linear  combination  is  used  with  two 
frequency  data  is  to  take  some  number  I  can  never  remember,  but  something  like  2.5  times  the  observable 
at  LI  minus  1.5  times  the  observable  at  L2.  It  is  very  possible  for  the  phases  of  the  multi-path  signals  at 
LI  and  L2  not  to  cancel,  but  to  add,  so  that  you  end  up  with  the  multi-path  effect  on  the  combined 
observable  being  as  much  as  4  times  what  it  was  on  either  of  the  single  frequencies.  So  I  offer  that  as  a 
caution. 

Another  thing  I  would  suggest,  as  I  have  done  many  times  in  the  past,  that  one  of  the  diagnostic  tools  that 
you  have  available  in  every  receiver  -  unfortunately,  it  has  been  poorly  supported  in  RINEX  -  is  to 
measure  the  amplitude  variability.  For  every  cycle  of  multi-path  you  go  through,  there  is  also  a  cycle  of 
amplitude.  And  in  voltage  units,  the  fractional  change  in  amplitude  is  precisely  the  same  thing  as  the 
phase  change  in  radians.  So  you  have,  in  a  sense,  an  orthogonal  measurement  that  you  can  use  as  a 
diagnostic  tool  for  identifying  multi-path  in  the  amplitudes. 

BILL  KLEPCZYNSKI:  Thank  you  very  much.  That  is  a  very  good  comment. 

KEN  JOHNSTON  (U.S.  Naval  Observatory):  1  have  one  comment,  Bill.  At  the  Observatory  we  will 
soon  be  getting  a  beam-forming  GPS  receiver,  a  16-element  receiver.  What  I  encourage  you  to  do  is  take 
that  receiver  and  compare  that  to  the  other  data  that  you  can  get  from  the  Naval  Observatory  and 
hopefully  report  on  that  next  year  at  the  PTTI  meeting. 
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